Could saturation effects be visible in a future electron-ion collider? 
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We expect to observe parton saturation in a future electron - ion collider. In this letter we discuss 
this expectation in more detail considering two different models which are in good agreement with the 
OO \ existing experimental data on nuclear structure functions. In particular, we study the predictions of 

saturation effects in electron - ion collisions at high energies, using a generalization for nuclear targets 
of the b-CGC model, which describes the ep HERA quite well. We estimate the total, longitudinal 
and charm structure functions in the dipole picture and compare them with the predictions obtained 
using collinear factorization and modern sets of nuclear parton distributions. Our results show that 
inclusive observables are not very useful in the search for saturation effects. In the small x region 
they are very difficult to disentangle from the predictions of the collinear approaches . This happens 
mainly because of the large uncertainties in the latter. On the other hand, our results indicate 
that the contribution of diffractive processes to the total cross section is about 20 % at large A and 
small Q 2 , allowing for a detailed study of diffractive observables. The study of diffractive processes 
becomes essential to observe parton saturation. 
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1^ ■ The search for signatures of parton saturation effects has been subject of an active research in the last years (for 
ly-j | recent reviews see, e.g. It has been observed that HERA data in the small x and low Q 2 region can be successfully 
described with saturation models Jl, [H, d, H, Q . Moreover, the experimentally measured total cross sections present 
the property of geometric scaling [7| and the observed supression of high px hadron yields at forward rapidities in dAu 
collisions at RHIC [8( follows the behavior predicted in the framework of the color glass condensate (CGC) formalism 
II1HI10- All these results provide some evidence for saturation at HERA and RHIC. However, more definite 
conclusions are not possible due to the small value of the saturation scale in the kincmatical range of HERA and 
■ due to the complexity of dAu collisions, where we need to consider the substructure of the projectile and the target, 
as well as the fragmentation of the produced partons. So far, other models (without saturation included) are able 
to describe the same set of data (see e.g. Refs. [3, [111). I n 01 "der to discriminate between these different models 
and test the CGC physics, it would be very important to consider an alternative search. To this purpose, the future 
electron- nucleus colliders are ideal 0, Il7l [To| , because they can probably determine whether parton distributions 
saturate or not. 

After the eRHIC was proposed [lTj , it became crucial to have some quan titative estimates of the impact of saturation 
effects on observables. Some of these estimates can be found in [H, M H EH- In particular, in [19( we calculated 
several inclusive observables {F^, F^, F^ A and their logarithmic slopes), while in [2(| the behavior of the diffractive 

structure function F^^ 1 and the diffractive cross section a dl ^ were studied in detail. An interesting conclusion from 
Ref. [2(| was related to the growth of R a = a dl ^ ja tot with the atomic number of the target, especially in the small 
x and low Q 2 region. This ratio could be as large as 0.3 — 0.4. This is very large compared to the correspondingratio 
in ep collisions which is of the order of 0.10 — 0.15. Such a large value of R a was anticipated long ago in Ref. [22[ but 
the investigation of this enhancement was not further developed. Literally, our result implies that in 30 to 40 % of 
the eA collisions the nucleus escapes intact! This phenomenon is so spectacular that it deserves further investigation. 

Our goal in this letter is twofold. Firstly, to improve our previous estimates for the inclusive observables at eA 
colliders considering a more realistic model for the nuclear dipole cross section which allows us to describe the scarce 
experimental data on nuclear structure functions. This can give us more confidence about the behavior predicted 
by the saturation model for F^, F^ and F£' A at the small x. Secondly, to compare these predictions with those 
obtained using collinear factorization and the parametrizations of the nuclear parton distributions. These improved 
calculations should allow us to give a partial answer to the question posed in the title of the paper. 
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In order to study the behavior of the observables of deep inelastic scattering (DIS) at small x and to include 
saturation effects, it is useful to describe the photon-hadron scattering in the dipole frame, in which most of the 
energy is carried by the hadron, while the photon has just enough energy to dissociate into a quark-antiquark pair 
before the scattering. This description contrasts with the usual description in the infinite momentum frame of the 
hadron, based on collinear factorization, where the photon scatters a sea quark, which is typically emitted by the 
small- x gluons in the hadron. The QCD description of DIS at small x can be interpreted as a two-step process p3l |. 
The virtual photon (emitted by the incident electron) splits into a qq dipole, with transverse separation r, which 
subsequently interacts with the target. In terms of virtual photon-target cross sections <jt,l for the transversely and 
longitudinally polarized photons, the nuclear F 2 structure function is given by F A (x, Q 2 ) = jjr^, — Ctou with [23| 

otot = o T + cr L and otx = J d 2 r dz \^t,l{v, z,Q 2 )\ 2 a A ip (x,r), (1) 

where ^t,l is the light-cone wave function of the virtual photon and o A ip is the dipole nucleus cross section describing 
the interaction of the qq dipole with the nucleus target. In equation (fTJ) r is the transverse separation of the qq pair 
and z is the photon momentum fraction carried by the quark (for details see e.g. Ref. [Hj])- Moreover, the nuclear 
longitudinal structure function is given by F A (x,Q 2 ) — (Q 2 /47r 2 a em ) ol and the charm component of the nuclear 
structure function ' (x, Q 2 ) is obtained directly from Eq. |T]) isolating the charm flavor and considering m c = 1.5 
GeV. 

The main input for the calculations of inclusive observables in the dipole picture is <J A ip {x, r) which is determined 
by the QCD dynamics at small x. In the eikonal approximation it is given by: 

a A p {x,r)=2 J d 2 bN A {x,r,b) (2) 

where Af A (x, r, b) is the forward scattering amplitude for a dipole with size r and impact parameter b which encodes 
all the information about the hadronic scattering, and thus about the non-linear and quantum effects in the hadron 
wave function. It can be obtained by solving the BK (JIMWLK) evolution equation in the rapidity Y = ln(l/x) 
Many groups have studied the numerical solution of the BK equation [25] , but several improvements are still necessary 
before using the solution in the calculation of observables. In particular, one needs to include the next-to-leading 
order corrections into the evolution equation and perform a global analysis of all small x data. It is a program in 
progress (for recent results see [26| ) . In the meantime it is necessary to use phenomenological models for Af A which 
capture the most essential properties of the solution. 

In [l9|, [2(| we assumed that the impact parameter dependence of N A can be factorized as N A {x, r, b) = 
Af A (x,r)S(b) and proposed a generalization for the nuclear case of the IIM model @ which was, at that time, 
the most sophisticated one. In other words, we have disregarded the impact parameter dependence of the scattering 
amplitude and assumed that the dipole nucleus scattering^ amplitude is related to the dipole proton one by a simple 
modification in the saturation scale: Q 2 p (x) — > Q 2 A — As x Q 2 p (x). This naive model is useful to obtain some idea 
about the magnitude of the saturation effects in eA processes in comparison to the linear case. However, in order to 
get more reliable predictions we should use a phenomenological model which describes the current scarce experimental 
data on the nuclear structure function as well as includes the impact parameter dependence in the dipole nucleus 
cross section. A model which satisfies these requirements was proposed some years ago in Ref. [27j . In this model 
the forward dipole-nucleus amplitude was parametrized as follows 



M A (x,r,b) = 1 -exp 



\AT A (b)oV dip {x,r 2 ) 



(3) 



where TU(b) is the nuclear profile function, which is obtained from a 3-parameter Fermi distribution for the nuclear 
density normalized to unity (for details see, e.g., Ref. [IH). The above equation, based on the Glauber-Gribov 
formalism 29], sums up all the multiple elastic rescattering diagrams of the qq pair and is justified for large coherence 
length, where the transverse separation r of partons in the multiparton Fock state of the photon becomes a conserved 
quantity, i. e. the size of the pair r becomes eigenvalue of the scattering matrix. It is important to emphasize that for 
very small values of x, other diagrams beyond the multiple Pomeron exchange considered here should contribute {e.g. 
Pomeron loops) and a more general approach for the high density (saturation) regime must be considered. However, 
we believe that this approach allows us to estimate the magnitude of the high density effects in the RHIC and LHC 
kinematic range. 

In [13] the author has assumed that a v dip was given by the GBW model. However, in the last years an intense 
activity in the area resulted in more sophisticated dipole proton cross sections, which had more theoretical constraints 
and which were able to give a better description of the more recent HERA data [l(| EH E3, EH, US HH ■ In what 





FIG. 1: Ratio Rf 2 = *2F^ / AF® for A = Pb and A — Ca. Although joined by a line the results of the saturation models have 
been computed for the ((x), (Q 2 )) of the data points from the E665 Collaboration [3^ ]. 

follows we will use the b-CGC model proposed in Ref. [3(J, which improves the IIM model [|[ with the inclusion of 
the impact parameter dependence in the dipole proton cross sections. The parameters of this model were recently 
fitted to describe the current HERA data [3lj]. Following [3(| we have: 
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with Y = ln(l/a;) and k = x"(7s)/x'(7s)j where \ is the LO BFKL characteristic function. The coefficients a and b 
are determined uniquely from the condition that Af p (x,r), and its derivative with respect to rQ Sl are continuous at 
rQ s = 2. In this model, the proton saturation scale Q s _ p now depends on the impact parameter: 



V X 



exp 



25, 
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(0) 



The parameter i?ccc was adjusted to give a good description of the ^dependence of exclusive J/ip photoproduction. 
Moreover the factors A/o and 7 S were taken to be free. In this way a very good description of F2 data was obtained. 
The parameter set which is going to be used here is the one presented in the second line of Table II of [3l|: 7^ = 0.46, 
B CGC = 7.5 GeV~ 2 , Af a = 0.558, x = 1.84 x lO" 6 and A = 0.119. 

In Fig. Q]we compare with the E665 data the predictions for the ratio Rp 2 = 2F^/AF® obtained using the 
b-CGC model as input in our calculations. For comparison we also show the predictions obtained using the GBW 
model. As expected, both models fail to describe the large x region. However, both models describe quite well the 
scarce experimental data in the region of small values of x and low Q 2 . The basic difference between the predictions 
occur in the normalization of the ratio, with b-CGC model predicting a large nuclear effect at small x. It is important 
to emphasize that after the choice of the model for the dipole proton cross section our predictions are parameter free. 
Since our model for the dipole nuclear dipole cross section, Eq. ([3]), describes reasonably well the experimental data 
for the nuclear structure function, we feel confident to compare its predictions with those obtained using the collinear 
factorization as well as in extending it to the calculation of the longitudinal and charm structure functions. 

In Fig. [5] we present our predictions for the x dependence of the ratio Rf 2 for two different values of the photon 
virtuality Q 2 . The two saturation models predict a similar behavior in the small x region, yielding values of w 0.7 
for x = 10~ 4 and Q 2 — 2.5 GeV 2 . This value grows to sa 0.8 for Q 2 = 10 GeV 2 For comparison we also show 
the predictions obtained using the collinear factorization and nuclear parton distributions resulting from the global 
analysis of the nuclear experimental data using the DGLAP evolution equation (for a recent discussion see [33[). Here 
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FIG. 2: Comparison between the predictions for Rf 2 of the saturation models (b-CGC and GBW) and the collinear ones (DS 
and EPS08) for two different values of Q 2 and A = Ph. 

we consider two different sets of nuclear parton distributions: the DS [34| and EPS 08 [35| nuclear parametrizations. 
They represent a lower and an upper bound for the magnitude of the nuclear effects, respectively. In particular, in the 
analysis of Ref. [35[ the authors include data on hadron production in the forward region at RHIC ||. As a result the 
amount of nuclear shadowing is considerably larger than the one obtained in other parametrizations. It is important 
to emphasize that the collinear predictions do not include dynamical effects associated to non-linear (saturation) 
physics, since they are based on the linear DGLAP dynamics. Therefore, the comparison between the saturation 
and collinear predictions could, in principle, tell us if the observable considered can be used to discriminate between 
linear and non-linear dynamics. The results shown in Fig. [5] demonstrate that this is not the case of the nuclear 
structure function. Although the predictions of the two collinear models are similar for large x (> 10 -2 ), where there 
exist experimental data, the predicted behavior at small x is very distinct. This difference is a consequence of the 
choice of different experimental data used in the global analysis and of the assumptions related to the behavior of 
the nuclear gluon density [HI]. The difference between the collinear predictions is so large at small x that it is not 
possible to extract any information about the presence or not of new QCD dynamical effects from the study of the 
nuclear structure function. 

Let us now consider the behavior of the nuclear longitudinal structure function. It is believed that this observable 
can be used to constrain the QCD dynamics at small x in ep collisions at HERA (See, e.g. [36j]). One of our goals is 
to verify if this assumption is also valid in eA collisions. Our results are show in Fig. [3] The predictions of saturation 
models are, as in the case of Rf 2 , contained within the uncertainty range of the collinear models. On the other hand, 
the difference between the two collinear models is bigger. In this case the results for large x are distinct too. This 
is due to the strong dependence of on the nuclear gluon distribution, which is very different in the two models 
considered (33j . 

Finally, we show in Fig. 2]our results for the ratio Rc = 2F^ A j AF2 ,D , which is determined by the charm component 
of the nuclear structure function. In this case the saturation model predictions and the DS collinear ones are similar, 
while the EPS08 one is very distinct. Similarly to the Rl case, this behavior is associated to the large magnitude 
of nuclear effects present in the nuclear gluon distribution predicted by the EPS08 parametrization, which has the 
strongest shadowing. We did not include in the plots the somewhat older but very well known EKS parametrization. 
Its predictions for Rp 2 , Rl and Rc would always lie between the GBW and EPS08 curves. So we can say that, even 
neglecting the EPS08 parametrization (which is the most extreme), the predictions of saturation models overlap with 
those from collinear models. 

Let us summaryze our results for inclusive observables F<f, F^ and F% , which should be measured in the first 
run of a future electron-ion collider. The conclusion that we can draw from the previous figures is that the two dipole 
models yield similar predictions for the structure functions and they fall inside the range of predictions of the EPS 
and DS parametrizations. The latter differ among each other because of the large freedom inherent to global data 
analysis. 

Although these observables are sensitive to saturation effects, as shown in [lj| by the comparison between the 
predictions of the full (linear + non-linear contributions) dipole models and their corresponding purely linear versions, 
it is not yet possible to draw any firm conclusion concerning the QCD dynamics from inclusive quantities. This is a 
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FIG. 4: Comparison between the predictions for Rc = ^F^' /AF^' of the saturation models (b-CGC and GBW) and the 
collinear one (DS and EPS08) for two different values of Q 2 and A = Pb. 



pessimistic partial answer for the question posed on the title, which implies that in order to discriminate the saturation 
effects wc should consider less inclusive observables. 

An alternative can be the study of the logarithmic slopes of the structure functions, as already discussed in fiol. f3Tj] . 
As pointed out in 37], the logarithmic slope of can be useful to address the boundary between the linear and 
saturation regimes, due to the A dependence present in the nuclear saturation scale (See also [lj|). Another possibility 
is the study of observables measured in diffractive deep inelastic scattering (DDIS), since the total diffractive cross 
section is much more sensitive to large-size dipoles than the inclusive one [2LI20I]. Basically, the saturation effects screen 
large-size dipole (soft) contributions, so that a fairly large fraction of the cross section is hard and hence eligible for 
a perturbative treatment. This was the main motivation of Ref. (2fj| . where we have computed observable quantities 
like R a — <J di ff/atot an d F^ 3 ^ in the dipole picture using a naive saturation model. One of the main conclusions 
was related to the growth of R a with the atomic number of the target, especially in the small x and low Q 2 region, 
where the ratio could be as large as 0.3 — 0.4. Although we postpone for a future publication a detailed analysis of 
the diffractive observables using the b-GCC model discussed in this letter, it is important to verify if the diffractive 
contribution for the total cross section is still large when estimated using a more realistic saturation model. Following 
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FIG. 5: Ratio of diffractive to total cross sections, R a , as a function of x with the b-CGC and GBW models. Long dashed 
and solid lines are for Pb targets at Q 2 = 5.0 and 1.0 GeV 2 , respectively. Dot-dot-dash and dot-dash lines are the same for Ca 
targets. 



[181 ] we estimate the total diffractive cross section as follows 

<r#ff=<rE + °? and a£ r = \ J d 2 r dz \* T , L {r, z, Q 2 )\ 2 J d 2 b {^\ (7) 

where da A ip /d 2 b — 2Af A (x, r, b), with J\f A given by Eq. ([3]). In Fig. [5] we present our predictions for the ratio R a as 
a function of x for two values of Q 2 and the atomic number considering the b-CGC and GBW models. In comparison 
with our previous results (2p| . we can observe that the inclusion of the impact parameter in the saturation model 
reduces the ratio in approximately 30 %. However, the values predicted are still large, especially for small Q 2 and 
large A. Our new predictions imply that the ratio for ePb collisions is a factor two larger than observed in ep collisions 
at HERA. Consequently, the study of diffractive processes should be an easy task in an eA collider, allowing for a 
detailed analysis of the QCD dynamics. Since the study of inclusive observables is not very illuminating, the study 
of diffractive processes becomes fundamental in order to constrain the QCD dynamics at high energies. 

As a summary, in this letter we have studied the predictions of saturation physics for electron - ion collisions at 
high energies, using a generalization for nuclear targets of the b-CGC model which describes the ep HERA quite well. 
We have estimated the nuclear structure function F A (x, Q 2 ), as well as the longitudinal and charm contributions and 
compared with the predictions obtained using collinear factorization and distinct sets of nuclear parton distributions. 
The basic idea is to compare the predictions from non-linear and linear QCD dynamics and verify if the experimental 
analysis of these observables in the future electron - ion collider could reveal the presence of saturation physics, as well 
as constrain the behavior of the saturation scale. Our results indicate that the inclusive observables are not adequate 
for this purpose due to the large uncertainty present in the collinear predictions at small x. On the other hand, it 
is expected that the collinear factorization formalism fails to describe diffractive eA processes, while the saturation 
formalism remains valid. The study of diffractive processes becomes essential to observe the saturation physics. In 
this letter we have estimated the contribution of these processes for the total cross section and demonstrated that it 
is about 20 % at large A and small Q 2 , allowing for a detailed study of diffractive observables. Our results motivate 
a more detailed study of the diffractive structure function as well as the leptoproduction of vector mesons in eA 
processes. These studies are fundamental in order to get a definitive answer for the question posed on the title. 
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